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Abstract  The  characteristics  of  oxy-coal  combustion  for  a  swirl  burner  with  a  specially  designed  preheating 
chamber  are  studied  numerically.  In  order  to  increase  the  accuracy  in  the  prediction  of  flame  temperature  and  igni¬ 
tion  position,  eddy  dissipation  concept  (EDC)  model  with  a  skeletal  chemical  reaction  mechanism  was  adopted  to 
describe  the  combustion  of  volatile  matter.  Simulation  was  conducted  under  six  oxidant  stream  conditions  with  dif¬ 
ferent  02/N2/C02  molar  ratios:  21/79/0,  30/70/0,  50/50/0,  21/0/79,  30/0/70  and  50/0/50.  Results  showed  that  02  en¬ 
richment  in  the  primary  oxidant  stream  is  in  favor  of  combustion  stabilization,  acceleration  of  ignition  and  increase 
of  maximum  flame  temperature,  while  the  full  substitution  of  N2  by  C02  in  the  oxidant  stream  delays  ignition  and 
decreases  the  maximum  flame  temperature.  However,  the  overall  flow  field  and  flame  shapes  in  these  cases  are  very 
similar  at  the  same  flow  rate  of  the  primary  oxidant  stream.  Combustion  characteristics  of  the  air-coal  is  similar  to 
that  of  the  oxy-coal  with  30%  O?  and  70%  C02  in  the  oxidant  stream,  indicating  that  the  rear  condition  is  suitable 
for  retrofitting  an  air-coal  fired  boiler  to  an  oxy-coal  one.  The  swirl  burner  with  a  specially  designed  preheating 
chamber  can  increase  flame  temperature,  accelerate  ignition  and  enhance  burning  intensity  of  pulverized  coal  under 
oxy-coal  combustion.  Also,  qualitative  experimental  validation  indicated  the  burner  can  reduce  the  overall  NOv 
emission  under  certain  02  enrichment  and  oxy-coal  combustion  conditions  against  the  air-coal  combustion. 
Keywords  oxy-coal,  eddy  dissipation  concept  model,  coal,  swirl  burner 


1  INTRODUCTION 

Conventional  pulverized  coal  boilers  use  air  as 
oxidant  and  thus  the  flue  gas  is  mainly  consisted  of 
nitrogen  (N2)  and  carbon  dioxide  (C02)  while  C02  is 
less  than  20%  in  volume.  The  capture  of  C02  from 
such  dilute  mixtures  using  amine  stripping  or  other 
techniques  is  expensive.  In  oxy-coal  combustion, 
oxygen  (02),  typically  of  greater  than  95%  purity,  is 
used  as  the  oxidant,  and  C02  concentration  in  the  flue 
gas  can  be  over  90%.  However,  for  flame  temperature 


control,  C02  dilution  into  the  oxidizer  stream  is  nec¬ 
essary  [1,2].  Such  a  manner  could  delay  the  ignition 
and  affect  combustion  process  of  pulverized  coals.  Thus, 
burners  with  high  performance  in  ignition  and  com¬ 
bustion  are  still  demanded  for  oxy-coal  combustion. 

In  our  previous  studies  [3,  4],  a  high  temperature 
air  combustion  (HTAC)  burner  was  designed.  The 
burner,  as  shown  in  Fig.  1,  is  specially  equipped  with 
a  preheating  chamber  and  the  pulverized  coal-containing 
primary  air  stream  enters  the  chamber  from  an  off-centre 
duct.  With  such  a  structure,  internal  high  temperature 
flue  gas  recirculation  can  be  established  in  the  chamber. 


Figure  1  Schematic  of  the  swirl  type  HTAC  burner  with  a  special  preheating  chamber 
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Experiments  and  applications  showed  that  the  HTAC 
burner  possesses  excellent  performance  in  combustion 
stabilization,  low  NO*  emission,  and  is  capable  to 
bum  low  volatile  content  (LVC)  coals  [3, 4] 

In  this  paper,  oxy-coal  combustion  with  the 
HTAC  burner  at  different  O2  concentrations  was  nu¬ 
merically  studied  and  eddy  dissipation  concept  (EDC) 
model  with  a  skeletal  chemical  reaction  mechanism 
was  adopted  to  describe  the  combustion  of  volatile 
matter.  The  results  were  compared  with  that  of  air 
combustion.  Simulation  was  also  conducted  for  the 
CE-enriched  air  combustion  in  which  N2  in  the  air 
stream  was  partially  replaced  by  O2  such  that  the  O2 
concentration  in  the  oxidant  stream  was  higher  than 
21%.  The  maximum  flame  temperature  and  flame 
temperature  contour  were  examined  to  assess  the  fea¬ 
sibility  of  boiler  retrofit  from  air-fuel  combustion  to 
oxy-fuel  combustion  with  the  minimal  cost. 

2  NUMERICAL  SIMULATION 

The  numerical  simulation  was  based  on  an  Eule- 
rian-Lagrangian  formulation  of  the  3D  gas-solid  two 
phase  flow,  using  the  commercial  CFD  code  FLEUNT 
6.3.  The  gas  phase  conservation  equations  were  solved 
in  Eulerian  frame  and  solid  phase  were  solved  in  La- 
grangian  frame,  and  the  coupling  between  two  phases 
was  done  through  interactive  source  terms.  Finite 
volume  pressure-based  solver  with  implicit  lineariza¬ 
tion  was  used.  The  SIMPLE  algorithm  was  used  for 
pressure-velocity  coupling.  The  space  derivatives  of 
the  diffusion  terms  were  discretized  by  the  central 
differencing  scheme,  and  the  stiff  non-linear  terms  were 
discretized  by  the  first  order  scheme. 

The  sub-models  for  turbulent  flow,  radiation,  par¬ 
ticle  tracking,  devolatilization,  volatile  matter  combus¬ 
tion,  char  combustion  are  introduced  as  follows. 

2.1  Model  for  turbulent  flow 

Several  two-equation  eddy  viscosity  turbulence 
models  and  Reynolds  stress  model  (RSM)  were  tested 
during  the  model  selection  for  the  combustion  of  pul¬ 
verized  coal.  The  comparison  between  simulation  results 
and  experimental  data  showed  RSM  has  good  per¬ 
formance  in  the  predictive  precision  and  stability  of 
convergence,  though  with  a  slightly  more  CPU  cost  [5]. 
Thus,  RSM  was  used  for  the  turbulent  flow,  and  stan¬ 
dard  wall-function  is  applied  to  resolve  near-wall  region. 

2.2  Model  for  radiation 

There  are  two  most  commonly  used  radiation 
models  for  particle  radiation,  namely  PI  model  [6,  7] 
and  discrete  ordinates  (DO)  model  [8,  9].  Due  to  its 
better  accuracy,  DO  method  was  chosen  in  spite  of 
higher  computing  cost.  Weighted  sum  of  gray  gases 
model  [10]  was  chosen  to  calculate  the  fluid  absorption 
coefficients. 


2.3  Model  for  particle  tracking 

The  status  and  position  of  coal  particles  were 
tracked  in  the  Lagrangian  frame.  The  size  of  coal  par¬ 
ticles  was  assumed  to  satisfy  the  Rosin-Rammler  dis¬ 
tribution.  Based  on  the  distribution,  the  coal  particle 
size  was  divided  into  10  intervals,  in  which  the  size  of 
the  particles  was  represented  by  their  averaged  value. 
Then,  in  each  cell  on  the  injection  plane,  the  particles 
with  the  same  size  were  evenly  divided  in  number  into 
3  groups.  The  totally  number  of  groups  of  coal  particles 
were  -12000.  In  each  group,  the  particles  were  re¬ 
garded  to  have  same  velocity,  temperature  and  other 
status,  behaving  the  same  way.  The  composition  and 
burning  rate  of  particles  of  each  group  change,  but  its 
size  remains  at  the  initial  one.  Such  a  treatment  greatly 
reduced  the  computational  cost  and  improved  the 
convergence.  The  trajectory  of  a  discrete  phase  parti¬ 
cle  was  predicted  by  integrating  the  forces  subjecting 
to  the  particles.  No  interaction  was  considered  between 
the  particle  groups,  while  the  interactions  between  the 
discrete  phase  with  the  continuous  phase  through  mass, 
momentum,  energy  and  species  sources  were  consid¬ 
ered.  In  the  simulation,  the  equations  for  continuous 
phase  and  discrete  phase  were  alternatively  solved  until 
solutions  of  both  phases  converge.  To  ensure  the  con¬ 
vergence  for  the  continuous  phase,  200-300  continuous 
phase  iterations  were  done  per  discrete  phase  iteration. 
The  dispersion  of  particles  due  to  turbulence  in  the  fluid 
phase  was  predicted  by  the  stochastic  tracking  model. 

2.4  Model  for  volatile  matter  combustion 

In  order  to  increase  the  accuracy  in  predicting 
flame  temperature  and  ignition  position,  the  EDC  model 
with  a  skeletal  chemical  reaction  mechanism  was  ap¬ 
plied  to  describe  the  combustion  of  volatile  matter. 
The  skeletal  mechanism  included  16  species  and  41 
reactions,  which  was  used  to  model  CH4/air  turbulent 
combustion  by  Yang  and  Pope  [11].  This  model  could 
avoid  the  over-prediction  of  temperature  and  prema¬ 
ture  ignition,  which  were  caused  by  fast  reaction  as¬ 
sumption  and  missing  of  intermediate  species  in  eddy 
dissipation  (ED)  model  and  mixture  fraction/PDF 
(possibility  density  function)  model.  In  the  framework 
of  FLUENT,  an  UDF  (user  defined  function)  code  was 
compiled  to  implement  the  above  features  used  in 
EDC  model,  based  on  the  suggestion  from  Byggstoyl 
and  Magnussen  [12]  and  Gran  et  al  [13]. 

The  EDC  model  was  developed  by  Magnussen 
and  Hjertager  [14]  in  1976  to  include  chemical  mecha¬ 
nisms  in  turbulent  reacting  flows,  based  on  eddy  dis¬ 
sipation  turbulence  energy  cascade  model  [15],  some¬ 
times  called  simply  eddy  cascade  model.  The  cascade 
model  connects  the  viscous  fine  structures  at  which 
combustion  takes  place,  and  the  larger  transporting 
eddies,  and  expresses  fine-structure  quantities  in  terms 
of  turbulence  energy  and  dissipation.  The  scales  of 
turbulence  are  assumed  to  be  continuously  distributed 
over  a  wide  time  and  space  spectrum.  Mechanical  energy 
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is  transferred  from  the  main  flow  to  large  eddies  and 
then  further  to  smaller  eddies.  The  larger  eddies  carry 
the  most  kinetic  energy.  Smaller  eddies  whirl  faster  but 
carry  less  energy.  The  smallest  eddies  have  the  highest 
frequency  and  the  largest  viscous  stresses.  Viscous 
friction  transfers  mechanical  energy  to  heat  which  is 
then  dissipated.  The  smallest  eddy  has  a  length  scale 
called  Kolmogorov  length  scale.  If  the  eddy  size  is 
smaller  than  this  scale,  no  turbulent  structure  exists. 
Based  on  this  idea,  the  concept  of  the  EDC  was  de¬ 
veloped:  There  is  a  region  in  which  the  reactants  can 
be  regarded  as  ideally  mixed  and  thus  chemical  reac¬ 
tion  kinetics  determines  the  rate  of  the  process,  while 
outside  this  region  the  reactants  are  not  mixed  and  do 
not  react.  Therefore,  each  computing  cell  is  separated 
into  two  regions:  reaction  region  and  inert  region,  and 
separately  called  as  fine-scale  structure  and  surround¬ 
ings.  Fine-scale  structure  is  the  smallest  eddy,  whose 
dimensions  are  small  in  one  or  two  directions,  however 
not  in  all  three.  These  fine  structures  have  the  charac¬ 
teristic  dimension  in  the  same  magnitude  as  Kolmogorov 
microscale,  and  mainly  located  in  the  highly  strained 
regions  between  larger  energy-rich  eddies  [16]. 

EDC  model  has  been  widely  used  in  gas  combus¬ 
tion  simulation  [17-19]  but  was  not  reported  for  coal 
combustion  with  detailed,  skeletal  or  reduced  chemistry 
for  volatile  matter  (VM)  combustion.  In  the  present 
study,  EDC  model  was  integrated  with  other  models  in 
FLUENT  software  with  a  skeletal  chemistry  for  VM 
combustion  and  it  is  expected  that  it  can  increase  the 
accuracy  in  prediction  of  flame  temperature  and  igni¬ 
tion  position. 

When  EDC  model  with  the  skeletal  chemical  re¬ 
action  mechanism  is  applied,  the  number  of  iteration 
of  direct  integration  is  huge.  To  reduce  the  CPU  cost, 
an  in-situ  adaptive  tabulation  (ISAT)  procedure  that 
was  proved  to  have  100-1000  times  speed-up  in  PDF 
modeling  of  a  turbulent  reacting  flow  by  Pope  [20] 
was  adopted  in  this  study. 

2.5  Definition  of  volatile  matter 

Although  great  efforts  have  been  paid  by  many 
researchers,  so  far  there  is  no  a  clear,  convictive  and 
popularly  accepted  way  to  define  the  volatile  matter 
(VM)  as  it  is  difficult  to  precisely  mimic  the  practical 
combustion  environment  and  measure  the  species  of 
volatile  matter  in  experiments.  According  to  the  sug¬ 
gestions  in  [21-23],  the  light  gases  CH4,  EE  and  CO 
were  regarded  as  the  main  VM  species.  Their  molar 
ratios  for  the  coal  used  were  defined  as  35  i  43  •  22, 
based  on  the  criterion  that  the  elemental  mass  should 
be  consistent  with  the  ultimate  analyses  and  total  heat 
released  should  be  the  same  as  the  higher  heating 
value  from  the  proximate  analysis. 

2.6  Model  for  devolatilization 


Currently,  there  are  several  commonly  used 


models  to  describe  volatile  matter  releases,  including 
the  first-order  reaction  model,  two-competing  model 
and  macromolecular  network  models  [24-31].  Since 
the  main  target  of  the  paper  is  to  study  oxy-coal  com¬ 
bustion  characteristic  with  the  HTAC  burner  at  differ¬ 
ent  O2  concentrations,  the  rather  complicated  macro- 
molecular  network  model  was  not  necessary  and  the 
two-competing  model  [31]  was  chosen.  The  model  can 
be  expressed  as  follows: 


-ma 


+  a2R2)  exp  -f  (tf,  +R2)dt 


d  t 


(1) 


where  R\  and  R2  are  competing  kinetic  rates  control¬ 
ling  the  devolatilization  over  different  temperature 
ranges;  mY(t )  is  volatile  yield  up  to  time  t,  kg;  is 
initial  particle  mass,  kg;  a\  and  a2  equal  to  0.3  and  1.0 
respectively;  ma  is  ash  content  in  the  particle,  kg;  fw,o 
is  the  initial  moisture  fraction  of  the  particle. 

The  two  kinetic  rates  are  expressed  in  the  Ar¬ 
rhenius  form  as 


Rx  =  A1  exp 
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According  to  Ref.  [31],  A\  and  E\  are  set  to 
2.0 xlO5  s_1  and  1.04 xlO8  J-mof1  respectively;  A2  and 
E2  are  set  to  1.3><107  s_1  and  1.67xl08  J-moE1  respec¬ 
tively  (built-in  parameters  in  Fluent  6.3). 


2.7  Model  for  char  surface  combustion 


The  kinetic/diffusion  surface  reaction  rate  model 
[32]  was  applied  to  describe  char  combustion,  assum¬ 
ing  that  the  surface  rate  of  particle  was  determined 
either  by  kinetics  of  surface  reaction  or  by  diffusion 
rate  of  oxygen.  The  model  can  be  represented  by 

dm  pRTJ^  1  (3) 

At  ^  Mw>ox  l/i?diff  +1/ Rc 

_o 

where  p  is  the  density  of  char,  kg-m  ;  Ap  is  the  area  of 
particle  surface,  m2;  R  is  the  universal  gas  constant, 
J-moE^K  ;  Tx  is  gas  temperature,  K;  Yox  is  the  oxi¬ 
dant  mass  fraction;  Mw?ox  is  the  molecular  mass  of  oxi¬ 
dant,  kg-kmol”1.  Rdiff  is  the  diffusion  reaction  rate  co¬ 
efficient  and  Rc  is  the  kinetic  or  chemical  reaction  rate 
coefficient.  Rdiff  and  Rc  are  expressed  as 
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Rc  =  C2  exp 
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where  dp  is  the  diameter  of  the  particle,  m;  Tp  is  the 
temperature  of  particle,  K;  Tg  is  the  surrounding  gas 
temperature  K;  Ea  is  the  reaction  activation  energy, 
J-kg^-moE1,  Ci  is  mass  diffusion-limited  rate  constant 
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(b)  Grid  of  vertical  section  the  burner  (c)  Grid  in  the  burner  front  section 

Figure  3  The  geometry  and  girds  of  burner 


and  C2  is  the  kinetics-limited  rate  pre-exponential 
factor.  The  constants  Cu  C2  and  Ea  are  5xl0”1-  kg-s”1, 
0.002  kg-s”1  and  7.9><lo"  Enrol”1  respectively  [32]. 

2.8  Structure  and  grid 

The  structure  and  dimensions  of  the  burner  and 
furnace  used  in  the  simulation  are  shown  in  Fig.  2. 
The  secondary  air  vane  angle  was  60°. 

The  hexahedral  type  grid  was  used  to  mesh  the 
computational  domain,  as  shown  in  Fig.  3.  This  type 
of  grid  not  only  was  more  stable  to  convergence,  but 
also  needed  much  less  mesh  number  than  other  types 
of  grids.  A  grid  independence  check  for  Case  1  was 
performed  in  order  to  properly  set  the  gird  number. 
First,  stable  solutions  were  obtained  under  the  same 
condition  but  with  different  grid  numbers.  Then,  a 
monitoring  point  was  selected  and  the  temperature 
difference  between  this  point  and  the  primary  air  inlet 
were  calculated.  When  the  temperature  difference  was 


nearly  independent  with  the  grid  number,  the  solution 
was  regarded  to  be  independent  of  grid  number.  In  the 
present  study,  a  monitoring  point  was  located  on  the 
centerline  of  the  furnace  and  50  cm  away  from  the 
burner  exit.  Shown  in  Fig.  4,  the  temperature  differ¬ 
ence  became  nearly  unchanged  (with  relative  error 


grid  number  ( N ) 

Figure  4  Grid  independence  check  for  Case  1 
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within  1%)  when  the  grid  number  was  larger  than 
35000,  indicating  that  further  increasing  the  grid 
number  would  not  increase  the  predictive  accuracy. 
Accordingly,  the  grid  number  of  40000  was  used  for 
Case  1 .  Due  to  the  similarity  in  numerical  simulation, 
this  gird  number  was  applied  to  other  cases  too. 

2.9  Coal  properties  and  operational  parameters 

The  proximate  and  ultimate  analyses  of  the  lean 
coal  used  in  simulation  are  listed  in  Table  1 .  The  main 
operational  parameters  of  the  burner  are  given  in 
Table  2.  The  primary  air  was  in  atmospheric  pressure, 
but  at  a  temperature  of  573  K.  The  excess  air  ratio  was 
remained  at  ~1.2  for  all  cases  and  correspondingly  the 
mass  flows  of  primary  and  secondary  air  were  ad¬ 
justed  at  different  O2  and  C02  concentrations.  In  the 


present  study,  six  cases  with  different  compositions 
and  mass  flow  rates  of  the  primary  and  secondary  airs 
as  listed  in  Table  3  were  simulated. 

3  RESULTS  AND  DISCUSSION 

3.1  Effect  of  02  concentration  in  the  primary  air 
on  overall  flow 

As  02  increase  in  the  primary  air  flow,  the  sec¬ 
ondary  air  flow  rate  and  thus  the  overall  flow  rate  of 
reacting  gas  decreases  in  some  extent,  as  exemplified 
by  Cases  1  to  3  in  Table  3.  It  is  then  interesting  to  check 
if  this  will  cause  significant  changes  in  the  overall 
flow  field.  Simulation  showed  that  the  velocity  fields 
of  six  cases  are  very  similar,  and  the  representative 
velocity  field  is  displayed  in  Fig.  5.  Four  reverse  flow 


Table  1  The  proximate  and  ultimate  analyses  of  coal  (air  dry) 

Proximate  analysis 

Ultimate  analysis/%  (by  mass) 

Mass  ratio/% 

Lower  heating  value 

C 

H 

N  S 

O 

Moisture 

Ash 

Volatile  Fixed  carbon 

/MJ-kg-1 

0.76 

41.90 

14.39  42.94 

17.70 

76.74 

4.14 

1.25  3.87 

14.00 

Table  2 

Operational  parameters  and  conditions  used  in  simulation 

Temperature/K 

Mass  flow  rate/kg  -s  1 

Pressure/kPa 

Mean  diameter/m 

Distribution 

primary  air 

512 

0.288 

101.3 

— 

— 

secondary  air 

593 

listed  in  Table  3 

101.3 

— 

— 

coal  particles 

512 

0.1347 

101.3 

5x1 0~6 

Rosin-Rammler 

Table  3  Composition  and  mass  flow  rate  for  the  primary  and  secondary  air 

Case 

^0/% 

<Pco2/% 

<pNJ% 

Fpri/kg-S  1 

V sec/kg '  S  ‘ 

1 

21 

0 

79 

0.288 

2.55 

2 

30 

0 

70 

0.288 

2.41 

3 

50 

0 

50 

0.288 

2.09 

4 

21 

79 

0 

0.288 

2.55 

5 

30 

70 

0 

0.288 

2.41 

6 

50 

50 

0 

0.288 

2.09 

Figure  5  Typical  map  of  velocity  vector  of  HTAC  burner  (Case  1) 
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Figure  6  Temperature  fields  with  different  02/N2/C02  composition  in  the  oxidant  stream 


zones  exist  in  the  area  adjacent  to  the  burner  including 
one  in  the  preheating  chamber,  one  outside  the  burner 
caused  by  swirl  secondary  air  and  two  side -wall  reverse 
flow  zones.  The  reverse  flow  in  the  preheating  cham¬ 
ber  entrains  hot  flue  gas  from  the  furnace  to  preheat 
the  primary  air,  and  thus  is  beneficial  to  flame  stabili¬ 
zation. 

3.2  Effect  of  primary  air  composition  on  tem¬ 
perature  field 

The  predicted  results  of  the  temperature  field 
with  six  O2/N2/CO2  compositions  in  the  oxidant  stream 
are  respectively  shown  in  Fig.  6.  Due  to  the  offset  of 
the  primary  gas  stream,  the  temperature  field  adjacent 
to  the  burner  is  asymmetric  and  in  most  cases,  flame 
center  locates  above  the  center  plane  of  the  burner.  Al¬ 
though  temperature  at  the  same  position  varies  with 
inlet  gas  composition,  the  flame  shapes  are  similar.  The 
results  indicate  that  difference  of  the  oxidant  stream 
composition  has  no  significant  effect  on  flame  shape, 
though  the  temperature  values  are  different.  This  is  an 
interesting  phenomenon. 

It  can  be  also  seen  the  temperature  field  for  the 
oxy-coal  combustion  with  30%  02  in  the  oxidant 
stream  is  close  to  the  one  in  air  condition.  Therefore, 
when  an  air-coal  coal  fired  boiler  is  retrofitted  to  an 
oxy-coal  boiler  for  carbon  capture  and  storage  (CCS) 
by  replacing  N2  by  C02,  increasing  02  concentration 
of  the  primary  air  to  30%  is  an  efficient  compensatory 
method,  and  this  is  consistent  with  the  literature  find¬ 
ings,  e.g.  in  Ref.  [2]. 

3.3  Effect  of  primary  air  composition  on  the 
flame  temperature 

The  maximum  flame  temperature,  denoted  as 
Tmax,  at  different  02,  N2  and  CO2  concentrations  in  the 
primary  oxidizer  stream  is  shown  in  Fig.  7.  It  can  be 


Figure  7  Maximum  flame  temperatures  with  different 
02/N2/C02  compositions  in  the  oxidant  stream 


seen  for  both  02-N2  and  O2-CO2  based  combustion, 
Tmax  increases  with  02  concentration,  while  at  the 
same  O2  concentration,  the  former  one  has  a  higher 
Tmax-  The  differential  temperature  between  02-N2  and 
02-C02  based  combustion  is  slightly  higher  at  a  lower 
O2  concentration  and  this  is  because  the  CO2  fraction 
is  correspondingly  higher  (as  shown  in  Table  3)  and 
thus  the  average  specific  heat  capacity  increases.  Con¬ 
sistent  with  results  shown  in  Fig.  4,  to  keep  Tmax  close 
to  the  one  in  air  based  combustion,  02  concentration 
in  the  oxidant  stream  should  be  increased  to  a  value  in 
the  range  of  30%~35%  when  oxy-coal  combustion  is 
adopted. 

The  increase  of  Tmax  in  the  02-enriehed  or 
oxy-coal  combustion  could  be  a  concern  for  the  addi¬ 
tional  NOx  formation.  On  one  hand,  higher  Tmax  does 
increases  thermal  NOx  formation.  However,  as  long  as 
the  average  temperature  in  the  main  combustion  zone 
is  less  than  1700  K,  thermal  NOx  occupies  only  a 
small  fraction  of  the  total  NOx  in  a  coal-fired  boiler 
[33].  On  the  other  hand,  fuel  NOx,  which  dominates 
the  total  NOx  formation  in  most  coal-air  combustion, 
could  be  remarkably  decreased  in  02-enriehed  com¬ 
bustion  when  the  ignition  of  the  coal  particles  is  kept 
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in  a  reducing  environment.  In  such  an  environment, 
02-enriehed  combustion  enhances  the  release  of  fuel-N 
into  volatile-N  in  high  temperature,  analogous  to  the 
low  NO*  combustion  feature  of  HTAC  [3, 4].  Thus,  the 
overall  NO*  formation  with  for  02-enriehed  or  oxy-coal 
combustion  is  not  necessarily  to  increases;  instead  it 
could  be  reduced  if  combustion  is  well  arranged.  The 
feasibility  of  overall  NO*  emission  reduction  is  vali¬ 
dated  by  experiments  introduced  in  the  later  section. 

3.4  Effect  of  primary  air  composition  on  ignition 
position 

So  far,  there  is  no  unified  definition  for  the  igni¬ 
tion  position  of  turbulent  coal  combustion.  In  the  pre¬ 
sent  study,  isothermal  curves  at  1173  K,  that  was  pro¬ 
posed  by  a  few  researchers  ( e.g .,  [34])  are  chosen  to 
define  the  ignition  position  of  pulverized  coal  in  the 
turbulent  jet  combustion.  The  numerically  defined 
ignition  positions  with  different  02  concentrations  in 
air-coal  and  oxy-coal  combustion  are  depicted  in  Fig.  8. 
With  the  increase  of  02  concentration,  ignition  occurs 
earlier  for  both  02-N2  and  02-C02  based  combustion. 
In  addition,  ignition  is  delayed  if  N2  replaced  by  C02 
in  the  oxidant  stream.  The  positions  at  which  ignition 
begins  (close  to  the  burner  exit)  for  the  cases  when 
(Po2  '■  <Pco.  =  30  :  70  and  cpQi  :  cpc0 ^  =  50  :  50  is 
close  to  that  for  the  air-coal  combustion,  while  the 
positions  at  which  ignition  begins  for  the  cases  when 
<p0  '■  cpco  =  21  -  79  is  much  delayed  than  that  for 
the"  air-coal  combustion.  The  results  are  mainly  attrib¬ 
uted  to  the  higher  specific  heat  capacity  of  C02  and  the 
lower  flame  temperature  in  the  oxy-coal  combustion. 


radial  distance  from  burner  centerline,  f?/mm 

Figure  8  The  ignition  position  under  different  02/N2/C02 
composition  in  the  oxidant  stream 

0  02-C02,  (Pq  =21%  ;  o  02-C02,  <Pq  =30%  ;  •  02-C02, 
(pQ  =50%;  ■  02-N2,  (pQ  =21%  ;  □  b2-N2,  (pQ  =30%; 
ffl  d2-N2,  (pQ  =  50% 

3.5  Qualitative  experimental  validation 

Some  experiments  on  the  characteristics  of 
02-enriehed  air-coal  combustion  were  conducted  in  a 
single  burner  test  installation  using  the  swirl  burner  as 
simulated.  The  rated  capacity  of  the  burner  was  about 
12  MW  in  thermal  capacity.  Several  coals  including 


the  one  used  in  above  simulation  were  tested,  includ¬ 
ing  a  lean  coal  with  volatile  (Fdaf=  14%)  and  heating 
value  ((2ar,net,P=  16800  kJ-kg”1)  similar  to  the  one  used 
in  simulation.  The  flow  rates  for  the  primary  air,  sec¬ 
ondary  air  and  over  fire  air  (OFA,  located  at  about  3.5 
m  downstream  from  the  burner)  were  about  0.65 
kg-s”1,  2.90  kg-s”1  and  0.68  kg-s”1  respectively.  The 
coal  feed  rate  was  about  0.42  kg-s”1.  The  overall  excess 
air  ratio  was  about  1.2.  All  air  streams  were  in  ambient 
temperature.  Because  of  the  experimental  limitations, 
the  furnace  preheated  at  around  700  K,  lower  than  that 
of  an  industrial  boiler  furnace.  During  the  experiments, 
in  order  to  save  02  consumption,  instead  of  increasing 
02  concentration  of  the  entire  primary  air  flow,  a  tube 
with  25  mm  inner  diameter  was  used  to  supply  an 
02-enriehed  02-N2  stream  into  the  preheating  chamber. 
The  flow  rate  of  the  02-enriehed  stream  is  constant  at 
9.5*  10”3  kg-s”1,  with  pure  02  stream  of  0-3.4*  10”3 
kg-s”1  balanced  by  air.  A  small  amount  of  assistant  oil 
(~20  kg-h”1)  was  also  used  since  flames  were  not  sta¬ 
ble  due  to  the  ambient  temperature  oxidizer  streams 
and  the  low  temperature  furnace.  NO*  measurements 
were  conducted  at  the  exhaust  section  of  the  flue  gas, 
after  combustion  remained  stable  for  about  5- 10  min. 

Though  the  experimental  conditions  were  not 
identical  to  the  ones  used  in  the  present  simulation, 
the  results  were  still  good  to  qualitatively  validate  the 
combustion  characteristics  and  NO*  emission  per¬ 
formance  of  the  swirl  burner.  Fig.  9  shows  that  NO* 
concentration  in  the  flue  gas  decreases  about  25% 
when  02  concentration  in  the  02-enriehed  stream 
increases  from  21%  to  50%.  Though  at  the  same  time, 
CO  concentration  increases  about  50%  due  to  the 
short  combustion  zone  in  the  experiments,  the 
results  are  sufficient  to  validate  the  hypothesis  that 
02-enriehed  or  oxy-coal  combustion  can  be  an  effec¬ 
tive  way  to  further  decrease  rather  than  increase  NO* 
formation  for  pulverized  coal  firing  when  the  proper 
burner  is  used.  The  reduction  mechanism  is  similar  to 
the  one  in  HTAC,  that  is  02-enriehed  or  oxy-coal 
combustion  can  be  in  favor  of  the  volatile  release 
before  coal  particles  disperse  in  the  oxidizing  envi¬ 
ronment  [4,  6]. 


02  concentration  in  the 
02-enriched  stream,  q>QJ% 


Figure  9  Experimental  data  on  variation  of  NO*  and  CO 
concentration  in  the  flue  gas  with  the  02  concentration  in 
the  02-enriched  flow 
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4  CONCLUSIONS 

Based  on  the  numerical  simulations  of  oxy-coal 
combustion  in  a  swirl  burner  with  a  specially  designed 
preheating  chamber,  it  was  found  that  there  is  obvious 
difference  between  oxy-coal  and  air-coal  combustion. 
The  substitution  of  N2  with  C02  in  the  same  volumet¬ 
ric  amount  will  decrease  the  maximum  flame  tem¬ 
perature,  and  delay  ignition  of  pulverized  coal  in  tur¬ 
bulent  stream,  but  the  overall  velocity  field  and  flame 
shapes  are  similar.  When  02  concentration  increases  to 
30%,  the  temperature  contour  adjacent  to  the  burner  is 
similar  to  the  one  in  air  combustion.  The  results  con¬ 
firmed  that  oxy-coal  combustion  with  30%  02  con¬ 
centration  is  feasible  for  retrofitting  the  existing 
burner  and  furnace  from  air-coal  combustion.  In  addi¬ 
tion,  eddy  dissipation  concept  model  with  skeletal 
chemical  reaction  mechanism  could  be  adopted  to 
describe  the  combustion  of  volatile  matter  in  pulver¬ 
ized  coal  combustion  simulation. 

The  swirl  burner  with  specially  designed  pre¬ 
heating  chamber  is  suitable  for  oxy-coal  combustion. 
The  structure  of  the  burner  is  beneficial  to  increase 
flame  temperature,  accelerate  ignition  of  pulverized 
coal  and  enhance  combustion.  The  simulation  results 
and  qualitative  experimental  validation  indicated  the 
burner  can  enhance  combustion  and  further  reduce  the 
overall  NOx  emission  under  certain  02-enriehed  and 
oxy-coal  combustion  conditions  against  the  air-coal 
combustion. 

NOMENCLATURE 

A  p  area  of  particle  surface,  m2 

A\  constant  (2. Ox  105  s_1) 

A 2  constant  (1.3xl07  s-1) 

Ci  mass  diffusion-limited  rate  constant  (5x1 0  12  kg-s% 

C2  kinetics-limited  rate  pre-exponential  factor  (0.002  kg-s  j 

dp  diameter  of  particle,  m 

Ea  char  surface  reaction  activation  energy,  J-moF1 

E i,  E2  activation  energy,  J  inol 

/w,o  initial  moisture  content  in  the  particle,  % 

MWtOX  molecular  mass  of  oxidant 

m  mass  of  char,  kg 

ma  ash  content  in  the  particle,  kg 

mVfi  initial  particle  mass,  kg 

mv(t)  volatile  yield  up  to  time  t,  kg 
R  universal  gas  constant,  J  •  mol- 1  -  K  1 

Rc  kinetic  or  chemical  reaction  rate,  kg-s  1 

diffusion  reaction  rate,  kg-s~‘ 

Ri,  R2  devolatilization  rate  over  different  temperature  range,  kg-s  1 

Tp,  Tg  particle  and  gas  temperature,  K 

IF  ambient  gas  temperature,  K 

V  mass  flow  rate,  kg-s 

F0x  oxidant  mass  fraction 

ai,a2  yield  factor 


p  density  of  char,  kg-m  3 

(p  volume  concentration,  % 
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